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The origin of the chemical elements in cluster cores 
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Metals play a fundamental role in ICM cooling processes in cluster cores through the emission of spectral lines. But when 
and how were these metals formed and distributed through the ICM? The X-ray band has the unique property of containing 
emission lines from all elements from carbon to zinc within the 0.1 — 10 keV band. Using XMM-Newton, the abundances 
of about 11 elements are studied, which contain valuable information about their origin. Most elements were formed 
in type la and core-collapse supemovae, which have very different chemical yields. Massive stars and AGB stars also 
contribute by providing most of the carbon and nitrogen in the ICM. Because feedback processes suppress star formation 
in the cluster centre, the element abundances allow us to directly probe the star formation history of the majority of stars 
that are thought to have formed between z = 2 — 3. The spatial distribution in the core and the evolution with redshift 
also provide information about how these elements are transported from the member galaxies to the ICM. I review the 
current progress in chemical enrichment studies of the ICM and give an outlook to the future opportunities provided by 
XMM-Newton's successors, like Astro-H. 
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1 Introduction 

After the 'Big Bang', the baryonic component of the Uni- 
verse mainly consisted of hydrogen and helium with traces 
of lithium and beryllium. The first metals with a higher 
atomic weight were produced in the first generation of stars, 
also referred to as Population III stars, which started the 
epoch of re-ionization. The nature of this stellar popula- 
tion is very uncertain, but it likely consisted of intermediate- 
mass and high-mass stars (Vangioni et al. 2011). Based on 
WMAP measurements, these stars formed around z 10, 
when the age of the Universe was about 500 Myr. Although 
these stars were the first to enrich the surrounding gas, the 
total contribution to the current day chemical composition is 
thought to be small, only 10"'' Zsoiar (Matteucci & Calura 
2005). 

The bulk of the enrichment probably occurred around 
z ^ 2 — 3 through the supernova explosions following 
major star bursts. A compilation of the universal star for- 
mation rate measurements as a function of redshift (Hop- 
kins & Beacom 2006) shows that the star formation rate 
peaked around z ^ 2 — 3 and declined slowly to low red- 
shift. However, these are averaged rates over both cluster 
and field galaxies. At a similar redshift, the Intra-Cluster 
Medium (ICM) in clusters of galaxies starts to form, which, 
together with feedback from Active Galactic Nuclei (AGN), 
quenched the star formation in these objects. Therefore, in 
the cluster of galaxies case, the star formation rate drops 
much faster than average (Gaboret al. 2010). In that respect, 
clusters of galaxies are a special environment, because their 
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enrichment is dominated by the products from the main star 
bursts at z ^ 2 — 3. 

Ferrara et al. (2005) have estimated that a major fraction 
(>90%) of the produced elements around z — 3 contributes 
to the enrichment of gas in a hot phase, instead of ending in 
newly formed stars. The ICM in clusters therefore contains 
the bulk of the metals produced in the cluster member galax- 
ies. In the period of the major star bursts, the main enrich- 
ment mechanisms that were transporting the metals from 
the galaxies to the surrounding medium were galactic winds 
driven by the supernova explosions and AGN uplifting of 
galactic gas. At a lower level, also other enrichment mech- 
anisms play a role. Metals can be ejected from the galaxies 
by galaxy-galaxy interactions, by sloshing motions of the 
hot ICM, and ram-pressure stripping of in-falling galaxies 
(see Schindler & Diaferio 2008 for a review). 



2 Origins of elements 

The bulk of the elements heavier than beryllium are pro- 
duced in the end phases of stellar evolution. The elements in 
the mass range between oxygen and silicon are mainly pro- 
duced by core-collapse supemovae, because the elements 
with masses higher than silicon that are produced in the core 
region of the massive star are compressed into the neutron 
star or black hole that is formed during the supernova event. 
Only the outer layers containing the lighter elements are 
ejected into the surrounding medium. The elements from 
silicon to nickel, however, are the main products of type 
la supemovae, which are exploding white dwarfs in binary 
systems. During such a supernova event, elements up to 
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Fig. 1 Expected abundances measured in a 120 ks XMM- 
Newton observation of Sersic 159-03 (bottom panel), which 
is a typical bright local cluster The statistical error bars 
were obtained from de Plaa et al. (2006). The estimates 
for the SNIa, SNcc, and AGB contributions are based on 
a sample of 22 clusters (de Plaa et al. 2007) and two ellip- 
tical galaxies (Grange et al. 201 1). The top panels show the 
typical range in SNIa and IMF models with respect to the 
statistical error bars in the observation. 

nickel are produced by explosive fusion and the white dwarf 
is disintegrated, which releases all the products in the sur- 
rounding medium. 

A few light elements, however, have a different origin. 
Carbon is thought to be mainly produced by massive stars 
during their lifetime and ejected through the stellar winds. 
Nitrogen is ejected by intermediate-mass stars in the Asymp- 
totic Giant Branch (AGB) phase. Since elements heavier 
than nickel are not produced by fusion, and are therefore 
less abundant, we do not discuss the origin of these metals 
in this paper 

Figure [T] shows the expected contributions of the most 
abundant metals that can be detected using the XMM-New- 
ton observatory. The bars indicate the relative contribution 
of each supernova type to the abundance of the elements. 
The fractions depend on the supernova models used. The 
effect of uncertainties in the supernova models and in the 
Initial-Mass Function (IMF) are indicated in the top panel. 

2.1 Type la supernovae 

Although type la supernovae are used as standard candles 
for cosmology, their progenitor and explosion mechanism 
are still poorly known. The common misconception is that 
a type la supernova occurs when a white dwarf goes over 
the Chandrasekhar mass. Instead, the explosive carbon fu- 
sion in the star is ignited just before it reaches this limit. In 
recent years, optical searches for supernovae have yielded 
hundreds of observations of type la's which revealed a sur- 
prising variety in their properties. It has proved to be very 
difficult to link the observed supernovae to possible progen- 
itors (see e.g. Howell 201 1 for a recent review). 



Currently, roughly three main progenitors are being con- 
sidered. The first is a 'classical' type la or also called 'Sin- 
gle Degenerate' (SD), where the companion of the white 
dwarf is a main sequence or red-giant star that accretes ma- 
terial on the white dwarf. The second is a scenario where 
two white dwarfs merge and the less massive white dwarf 
is accreted onto the more massive one, which is known as 
the 'Double Degenerate' (DD) scenario. And the last is a 
sub-Chandrasekhar channel, where a thick layer of helium 
builds up on the white dwarf's surface either by hydrogen 
burning or by a helium-rich donor. Because of the lower 
density of the last systems, they are thought to produce more 
intermediate-mass elements like silicon, sulfur, and calcium 
than during the deflagration in the typical SD and DD cases. 

Observations have shown that relatively luminous type 
la supernovae tend to occur in spiral galaxies, while the 
sub luminous are mainly found in elliptical galaxies with 
old stellar populations (e.g. Howell 2001; Sullivan et al. 
2006). This appears to indicate that the brighter 'prompt' 
supernovae explode early (<400 Myr) after the star burst, 
while the sub-luminous occur with a delay of a few Gyr. 
This can be explained in the DD scenario by the fact that 
it takes a longer time to evolve a star into a low-mass white 
dwarf compared to a high-mass white dwarf. The delay time 
explanation is also supported by high-redshift studies that 
show that at z = 1 type la supernovae are 12% more lu- 
minous and have less intermediate mass elements in their 
spectra than local SNIa's (e.g. Howell et al. 2007; Sullivan 
et al. 2009). However, also differences in initial metallicities 
or other unknown parameters could play a role. 

One of the few successes to combine theory and obser- 
vations in this field is the explanation of the observed Delay- 
Time Distribution (DTD) using stellar population synthesis 
studies. Supernova rate measurements in elliptical galaxies 
(and clusters of galaxies) show that the rate declines with 
~ t^^ after the initial star burst. Model calculations show 
the same level of decline for DD mergers, which is an indi- 
cation that this channel is dominating the SNIa rate at larger 
delay times (see Fig|2l Ruiter et al. 2009; Mennekens et al. 
2010). 

Due to the uncertainty in the nature of the SNIa progeni- 
tor and the poorly understood physics of the explosion (dur- 
ing the explosion, the burning front probably switches from 
being subsonic to supersonic), there is a large range of su- 
pernova la models predicting the metal yields. Well known 
examples are Iwamoto et al. (1999) and models by Bravo 
et al. (2012). The predicted yields for various elements can 
vary substantially for every model, which means that ac- 
curate abundance measurements can help to constrain the 
SNIa models. 

2.2 Core-collapse supernovae 

Although there is also a range of different models for core- 
collapse supernova yields in the literature (e.g. Woosley & 
Weaver, 1995; Tsujimoto et al., 1995; Chieffi & Limongi, 
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Fig. 2 Supernova type la delay time distributions for three 
different progenitor channels: white-dwarf mergers (DDS), 
Single-Degenerate Systems (SDS), and AM CVn systems 
as described in Ruiter et al. (2009). AM CVn systems are 
ultra-compact systems with a WD accretor and a He-rich 
donor. The top and bottom panel represent the result from 
two choices of initial parameters, but both are calculated 
for stellar populations in elliptical galaxies. The dotted line 
shows the observed decline of the type la supernova rate of 
~ where t is the time passed since the initial star burst. 
The decline is consistent with the DDS progenitor model. 



2004), the main parameters that determine the abundances 
measured are initial metallicity and the choice of the IMF. 
In order to determine the total contribution of core-collapse 
supernovae to the cluster enrichment, the model yields for 
each mass bin need to be weighted with the IMF in the fol- 
lowing way (Tsujimoto et al, 1995): 

„50 rT , \ 

M — ^ m 



^50 





where Mi{m) is the ith element mass produced in a main- 
sequence star of mass m. A standard Salpeter IMF corre- 
sponds to X = 1.35, but this equation can, of course, be 
altered to represent other IMFs. 

2.3 Intermediate-mass AGB stars 

Since intermediate-mass stars in their AGB phase are a ma- 
jor source of nitrogen, it is important to include them in the 
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Fig. 3 EPIC spectrum of the cluster 2A 0335H-096 with 
an exposure time of 130 ks. Adapted from Werner et al. 
(2006b). 

model as well. Since the yields depend on the initial mass 
of the main sequence star, we need to weigh the yields also 
with the IMF like in Equation[T] AGB star yields for differ- 
ent masses and initial metallicities are calculated by, for ex- 
ample, Karakas (2010). Together with the information from 
the core-collapse models, the AGB stars provide an addi- 
tional constraint to the initial metallicity and IMF of the 
stellar population. 

3 Abundance measurements in X-rays 

The soft X-ray band between 0.1 and 10 keV is very suit- 
able to measure abundances in hot plasmas, because it con- 
tains spectral lines of all elements from carbon to nickel. 
The hot ICM in clusters is particularly interesting, because 
the plasma is in collisional-ionization equilibrium, which 
simplifies the determination of the abundances. Moreover, 
the cluster ICM contains the integrated yield of billions of 
supernovae, providing a general picture of supernova yields, 
contrary to studies of a few individual supernovae through 
their remnants in our galaxy. The sensitivity of XMM-New- 
ton, and also Suzaku, allows us to measure abundances with 
enough accuracy to constrain supernova models. An exam- 
ple of a deep EPIC observation of the cluster 2A 0335H-096 
in Figure[3]shows the position and strength of the lines of the 
most abundant elements. The availability of these excellent 
data sets has triggered a number of successful abundance 
studies (see also the review by Werner et al. 2008). 

3.1 Enrichment in local clusters 

The first attempt to link measured abundances in the hot 
ICM in local clusters (z < 0.2) was performed using the 
ASCA satelHte (Mushotzky et al., 1996). Its instruments al- 
lowed the accurate detection of O, Ne, Mg, Si, S, Ar, Ca, Fe, 
and Ni for the first time. From this and later ASCA studies 



Copyright line will be provided by the publisher 



4 



J. de Plaa: The origin of tlie cliemical elements in cluster cores 



in 



Si 


s 


Ar 


Co 


1 

Fe 


Ni 








































SNIa 










. + .. 




SNcc 














SNIa + SNcc 

1 . . . 



15 20 25 

Atomic Number 



Fig. 4 Abundance ratios fitted with supernova yields from 
the WDD2 SNIa model (Iwamoto et al. 1999) and a SNcc 
model with an initial metallicity Z = 0.02 and a Salpeter 
IMF. The calcium abundance appears to be underestimated. 
It can not explain the Ar/Ca ratio measured in this XMM- 
Newton sample of 22 clusters. (Adapted from de Plaa et al. 
2007). 

(e.g. Fukazawa et al. 1998; Finoguenov et al. 2000) a gen- 
eral picture emerged of an ICM which was enriched early 
in the cluster evolution by core-collapse supernova and later 
by delayed type la supernovae. 

When more deep observations of clusters with XMM- 
Newton became available, de Plaa et al. (2007) performed 
an abundance study in a sample of 22 clusters with EPIC. 
With the obtained accuracy, it was possible to estimate the 
relative contributions of type la and core-collapse super- 
nova to the enrichment of the cluster by linearly combining 
the supernova yields for each type. For the most commonly 
used type la models (Iwamoto et al. 1999), this results in a 
SNIa contribution of about 30%, which is remarkably simi- 
lar to the SNIa/SNcc ratio of ^ 0.2 - 0.4 as observed in the 
optical band (e.g. Horiuchi et al. 2011) given the fact that 
de Plaa et al. (2007) ignored galactic evolution effects. Sur- 
prisingly, the calcium abundance did not fit the Iwamoto et 
al. (1999) models well (see Fig|4|i, while a model by Bravo 
et al. (1996), which also fitted to the Tycho supernova rem- 
nant (Badenes et al. 2006), resulted in a much better fit of 
the Ar/Ca ratio (de Plaa et al. 2007). This showed that X-ray 
spectra enable us to constrain type la supernova models. 

Several groups have performed similar studies with XMM- 
Newton and Suzaku data since then, but usually on fewer 
clusters and with varying numbers of elements. In a small 
sample, Sato et al. (2007) performed the same fit as de Plaa 
et al. (2007) using abundances measured with Suzaku XIS. 
Unfortunately, argon and calcium were not included in the 
fit, but their main conclusion was consistent with de Plaa 
et al. (2007) considering the uncertainty introduced by the 
limited number of measured elements. Recently, Bulbul et 
al. (2012) used a modified APEC model that is able to fit 



the supernova type la to core-collapse ratio directly to the 
spectra. They also report a type la contribution of 30^0%, 
which is consistent with previous work and optical data. It 
should be noted that this ratio depends highly on the as- 
sumed supernova models (De Grandi & Molendi, 2009), but 
it is reassuring that different groups find similar numbers in- 
dependently. 

3.2 The abundance of nitrogen and carbon 

In recent years, deep XMM-Newton observations of giant 
elliptical galaxies also allowed the detection of carbon and 
nitrogen. These two elements are not produced in large quan- 
tities in core-collapse supernovae, but in a variety of sources. 
The main origins of these elements are still subject of de- 
bate (see e.g. Romano et al. 2010), but likely candidates 
are metal-poor massive stars and intermediate-mass stars in 
their AGB phase. 

Since carbon and nitrogen have a relatively low atomic 
mass, their K-shell transitions are located in the soft X-ray 
band below 0.5 keV. It also causes the hnes to be stronger at 
lower temperatures < 1 keV (at temperatures of a few keV, 
a larger fraction of the carbon and nitrogen is fully ionised 
and not producing line emission). To detect the lines, it is 
best to observe cool elliptical galaxies with the Reflection 
Grating Spectrometer (RGS) aboard XMM-Newton. This 
instrument has sufficient effective area at low energies and 
a high spectral resolution to resolve the lines of these el- 
ements. Unfortunately, noise in the RGS CCD in the band 
where the carbon line is located prevents an accurate carbon 
measurement for weak sources. 

Werner et al. (2006a) and Grange et al. (2011) anal- 
ysed RGS data of the giant elliptical galaxies M87 (Werner), 
NGC 5044, and NGC 5813 (Grange) and detected nitrogen. 
They found that the nitrogen abundances in these objects are 
indeed much higher than expected based on supernova mod- 
els. From the high N/O ratio, it can be deduced that nitro- 
gen should still be produced in low- and intermediate-mass 
stars contrary to supernovae, because the star-formation rate 
in elliptical galaxies has declined very rapidly after the star 
bursts around 2^2 — 3. 

3.3 Intermediate redshift clusters 

Studying the chemical enrichment of clusters above z = 0.2 
becomes increasingly difficult, because of the lower flux 
of the sources. The Fe-K feature is usually the only fea- 
ture strong enough to resolve at higher redshifts. Therefore, 
studies at these large distances focus on the iron abundance 
evolution as a function of redshift. In a sample of 56 clus- 
ters with redshifts in the 0.3< z <1.3 range observed with 
XMM-Newton and Chandra, Balestra et al. (2007) found 
a significant decreasing trend of the iron abundance with 
redshift from roughly 0.4 solar at z = 0.3 to 0.2 solar at 
z = 1.3. It was, however, also clear that the scatter on 
the measured data points was relatively large. In a different 
sample of 39 clusters observed with XMM-Newton, Baldi 
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Fig. 5 Measured abundances versus redshift for a sample 
of intermediate redshift clusters. The top panel shows the 
abundances measured up to O.6R500 and the bottom panel 
the results when the core of the cluster is ignored. Adapted 
from Baldi et al. (2012). 

et al. (2012) recently reported that they could not confirm 
the trend found by Balestra et al. (2007), because of the 
high level of scatter in the data points (see Figure [Sj. Re- 
moving the cool cores in the spectral analysis just had a 
limited effect on the amount of scatter and did not allow 
them to draw firm conclusions. Since these measurements 
still only sample the tail of the SNIa delay-time distribu- 
tion that originates at z = 2 — 3, a strong trend may not 
be expected. Clearly, instruments with a much higher effec- 
tive area are needed to reduce the statistical errors on these 
measurements. 

4 Discussion and future prospects 

X-ray spectroscopy of clusters of galaxies has proved to be 
a useful tool to study the enrichment history of the hot ICM, 
but the field is still in its early stages of development. Using 
current instruments, it is possible to measure the abundance 
of about 10 elements with an accuracy of 20-30% if one 
includes systematic errors (De Grandi & Molendi, 2009). 
This appears to be a substantial systematic error, however, 
the uncertainty in the supernova yields, given by the spread 
in supernova models, is sometimes more than a factor of two 
for certain elements. This means that even with these sys- 
tematic uncertainties, the measurements can put constraints 
on supernova models and the IMF. 

There are two main sources of systematic errors in abun- 
dance measurements: the atomic data in spectroscopic codes 
and instrument calibration issues. A third one can be intro- 
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Fig. 6 Simulated Astro-H micro-calorimeter spectrum of 
Abell 1795 compared to an XMM-Newton EPIC MOS 
spectrum. The spectra were calculated using an updated ver- 
sion of the atomic database in SPEX including trace ele- 
ments, like CI, K, Ti, Cr, and Mn. 

duced when the multi-temperature structure in the gas is 
not taken properly into account in the spectral fit (Buote 
& Fabian 1998). Both the systematics in the atomic data 
and the calibration depend very much on the specific tran- 
sition and wavelength of the line. The silicon abundance, 
for example, is mainly determined by measuring the Si XIV 
lines at about 2.0 keV, which is very close to the Au-edge 
of the mirror and Si-edge of the detector This is a very 
challenging wavelength band to calibrate accurately. In this 
case, the Si XIV lines are well characterized in the atomic 
database, while there are also other lines from, for example, 
Fe XVII that have a much larger uncertainty in their oscil- 
lator strengths (de Plaa et al., 2012). It is very important to 
understand the systematics for each measured element indi- 
vidually to assess the quality of the supernova fits. 

In general, most of the measured abundances have a sys- 
tematic uncertainty that is lower than or not much higher 
than 20%, which is still much better than the scatter in the 
models. To beat down these systematic errors, investments 
have to be made in the development of spectral codes and 
their underlying atomic data. In addition, a high level of cal- 
ibration accuracy is to be pursued to improve the abundance 
accuracy both for current and future missions. 

4.1 Future prospects 

The next leap forward will be made using micro-calorimeter 
technology, which will allow spatially resolved high-reso- 
lution X-ray spectroscopy of extended sources. This type 
of sensor, which will be flown on Astro-H (Takahashi et 
al. 2010), has a typical resolution of a few eV across the 
soft X-ray band. With this technology, it will be possible 
to resolve lines from less abundant elements, like sodium, 
chromium, and manganese (see Fig |6|. Knowledge about 
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the abundances of these trace elements puts additional con- 
straints on supernova models. The Mn/Cr ratio, for example, 
appears to be very sensitive to the metallicity of the type la 
progenitor (Badenes et al. 2008) and the Na abundance is, 
like nitrogen, also correlated with the contribution of inter- 
mediate mass AGB stars and sensitive to the IMF. 

To measure trace elements like chromium and manga- 
nese in high resolution spectra, the spectral fitting codes 
that are used to fit the spectra need updates to their atomic 
data. Until recently, only the lines from elements with a high 
abundance (typically with even atomic numbers) where tak- 
en into account in the codes. The commonly used databases 
like ATOMDbQ (Smith et al. 2001) and the SPEX| Hne 
database (Kaastra et al. 1996) are being updated and are 
(partly) available in the spectral fitting codes. The system- 
atic errors in the atomic data, which can reach the 20% level 
in some important lines, remain an issue. Further invest- 
ments in laboratory measurements of lines and theoretical 
line calculations are needed to improve the accuracy. 

Until now, the measured abundances have been com- 
pared to individual supernova models. It would, however, 
be much more realistic to estimate the expected abundances 
from advanced binary population synthesis codes. Using these 
codes, it is possible to give the proper weight to the different 
type la progenitor scenarios and their respective yields. The 
high calcium abundance measured by de Plaa et al. (2007), 
for example, might be explained by sub-Chandrasekhar type 
la's, because explosive helium fusion produces more inter- 
mediate mass elements, like calcium. Predicted abundances 
using these detailed binary population synthesis can be di- 
rectly compared to the measured abundances in clusters. 
The high amount of measurable elements will substantially 
improve the quality of the test and therefore also our knowl- 
edge of chemical enrichment in clusters. 
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